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In this paper we study the gravitational field of a straight string generated from a class of nonlinear 
sigma models, specifically the Skyrme model without a twist and the Skyrme model with a twist 
(the twisted Skyrmion string). The twist term, mkz, is included to stabilize the vortex solution. 

To model the effects of gravity, we replace the Minkowski tensor, rf‘ v , in the standard Skyrme 
Lagrangian density with a space-time metric tensor, g ^, assumed to be static and cylindrically 
symmetric. The Einstein equations for the metric and field components are then derived. This work 
is still in progress. 


I. INTRODUCTION TO THE 0(3) NONLINEAR SIGMA MODEL 

A nonlinear sigma model is an A-component scalar field theory in which the fields are functions defining a mapping 
from the space-time to a target manifold^. By a nonlinear sigma model, we mean a field theory with the following 
properties^: 

(1) The fields, (j>(x), of the model are subject to nonlinear constraints at all points x € Afo, where Ado is the source 
(base) manifold, i.e. a spatial submanifold of the (2+1) or (3+l)-dimensional space-time manifold. 

(2) The constraints and the Lagrangian density are invariant under the action of a global (space-independent) 
symmetry group, G, on <j>(x). 

The Lagrangian density of a free (without potential) nonlinear sigma model on a Minkowski background space-time 
is defined to be£ 


£ = 2^2 T4b(</>) rf v dy.(j> A 9 v (j) B 


(1) 


where 7 ab{4> ) is the field metric, rj fJ,v is the Minkowski metric tensor, A is a scaling constant with dimensions of 
(length/energy) 1 / 2 and (f> = cp A is the collection of fields. Greek indices run from 0 to d— 1, where d is the dimension 
of the space-time, and upper-case Latin indices run from 1 to A. 

The simplest example of a nonlinear sigma model is the O(A) model, which consists of A real scalar fields, <j> A , 
with the Lagrangian density^ 


1 ^ d(p A d(j> A 

2A 2 Sab V dx» dx v 


where the scalar fields, (f> A , satisfy the constraint 


Sab 4 A t B = 1 


( 2 ) 


( 3 ) 


and Sab is the Kronecker delta. The Lagrangian density © is obviously invariant under the global (space independent) 
orthogonal transformations 0(A), i.e. the group of A-dimensional rotations^ 


b' A =oi <f B . 


( 4 ) 


One of the most interesting examples of a O(A) nonlinear sigma model, due to its topological properties, is the 
0(3) nonlinear sigma model in 1+1 dimensions, with the Lagrangian density^ 


£ = ^2 ^ d ^ 


( 5 ) 




2 


where p and v range over {0,1}, and <f> = (0 1 , </> 2 , </> 3 ), subject to the constraint </>•</>= 1, where the dot (.) denotes 
the standard inner product on R 3 . 

For a 0(3) nonlinear sigma model in any number d of space-time dimensions the target manifold is the unit sphere 
S 2 in R 3 , and p and v in the Lagrangian density © run from 0 to d— 1. A simple representation of <f> (in the general 
time-dependent case) is 


/sin/(f, r) sin 5 (f,r)\ 
<j) = sin/(i,r) cos g(t, r) 

\ cos f(t, r) ) 


( 6 ) 


where / and g are scalar functions on the background space-time, with Minkowski coordinates x^ = (t, r). In what 
follows, the space-time dimension d is taken to be 4, and so r is a 3-vector. 

If we substitute © into ©. the Lagrangian density becomes 

C = + f) ^ ^ 9 du9] (7) 
The Euler-Lagrange equations associated with in dTJ) are 

ff v d^d„f - (sin / cos /) tf v d^g d v g = 0 (8) 

and 

d^d v g + 2(cot /) d^f d v g = 0. (9) 


II. SOLITON SOLUTION 

Two solutions to the 0(3) field equations © and © are 
(i) a monopole solution, which has form 


(x/p\ 

cj) = r= \ y/p\ (10) 

W p ) 


where p = ( x 2 + y 2 + z 2 ) 1 / 2 is the spherical radius; and 
(ii) a vortex solution, which is found by imposing the 2-dimensional ” hedgehog” ansatz 

( sin/(r) sin< 7 (r)\ /sin f(r) sm(nO — x)\ 

sin/(r) cos g(r) = I sin/(r) cos {nO — \) (11) 

cos f{r) ) \ cos f(r) ) 

where r = ( x 2 + y 2 ) 1 ^ 2 , 0 = arctan (x/y), n is a positive integer, and \ is a constant phase factor. 

A vortex is a stable time-independent solution to a set of classical field equations that has finite energy in two spatial 
dimensions; it is a two-dimensional soliton. In three spatial dimensions, a vortex becomes a string, a classical solution 
with finite energy per unit length^. Solutions with finite energy, satisfying the appropriate boundary conditions, are 
candidate soliton solutions^. 

The boundary conditions that are normally imposed on the vortex solution (HID are /(0) = 7r and lim r _>.oo /(r) = 0, 
so that the vortex ’’unwinds” from <f> = —z to </> = z as r increases from 0 to oo. The function / in this case satisfies 
the field equation 


dr 2 dr 


— sin / cos / = 0 
r 


There is in fact a family of solutions to this equation satisfying the standard boundary conditions: 


sin / = 


2K 1 / 1 r n 
1 + I<r 2n 


( 12 ) 


(13) 
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or equivalently 


cos / = 


Kr 2n - 1 
Kr 2n + 1 


where K is positive constant. 

The energy density of a static (time-independent) field with Lagrangian density, £, J7J) is 


- £ = ~ ^2 d ^f d "f + ( sin2 /) ^ 


The energy density of the monopole solution is 

- C= 1 

A V 

while the energy density of the vortex solutions is 

AKn 2 r 2n ~ 2 
= A 2 (A> 2n + l) 2 ' 

In particular, the total energy 


E = 


(—£) dx dy dz, 


(14) 

(15) 

(16) 

(17) 

(18) 


of both the monopole solution and the vortex solutions is infinite. But the energy per unit length of the vortex 
solutions 

y = J J (-£) dx dy = (19) 

is finite, and does not depend on the value of K. (We use the same symbol for the energy per unit length and the 
mass per unit length, due to the equivalence of energy and mass embodied in the relation E = me 2 . Here, we choose 
units in which c = 1). 

This last fact means that the vortex solutions in the nonlinear sigma models have no preferred scale. A small value 
of K corresponds to a more extended vortex solution, and a larger value of K corresponds to a more compact vortex 
solution, as can be seen by plotting / (or —£) for different values of K and a fixed value of n. This means that the 
vortex solutions are what is called neutrally stable to changes in scale. As K changes, the scale of the vortex changes, 
but the mass per unit length, /j,, does not. Note that because of equation (|T9|) . there is a preferred winding number, 
n = 1, corresponding to the smallest possible positive value of y. 

Furthermore, it can be shown that the topological charge, T, of the vortex defined by 

T = £abc I f f +> A d x (j) B dy4> c dx dy (20) 

where £abc is the Levi-Civita symbol, is conserved, in the sense that <9 t T = 0 no matter what coordinate dependence 
is assumed for / and g in Q3D- So the topological charge is a constant, even when the vortex solutions are perturbed. 
Also, it is simply shown that for the vortex solutions 

T = -\ n I/(°o) - /(0)] = n (21) 

and so the winding number is just the topological charge. 

Because there is no natural size for the vortex solutions, we can attempt to stabilize them by adding a Skyrme term 
to the Lagrangian density^. 


III. SKYRMION VORTEX WITHOUT A TWIST: THE SKYRME TERM 

The original sigma model Lagrangian density (with the unit sphere as target manifold) is 

Cl = ^V2 T l tlU d 


( 22 ) 
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If a Skyrme term is added to (l22l) . the result is a modified Lagrangian density 


£2 = 7 tt "2 dn<t> . d u (j) - K s rj KX r]^" (d K (/) x d^) . ( d\(j) x d v tj>) 
ZA* v---✓ 


Skyrme term 


where the Skyrme term is the second term on the right hand side of 
With the choice of field representation equation (E^l) becomes 

£2 = Z— (rf v d^f d v f + sin 2 / if" d^g d v g) 


- K, 

The Euler-Lagrange equations generated by £ 2 , namely 

a, 

and 

d, 


2 sin 2 / ( if" d^f d v f) ( rf x d K g d x g ) - 2 sin 2 / {rT dj d v gf 


dC 2 \ 

dC 2 

d(d a f)J 

df 

dC 2 \ 

dC 2 

d(d a g)J 

dg 


= 0 


= 0 


(23) 


(24) 


(25) 


(26) 


will be considered more closely in Section IV. The mass per unit length of a static vortex solution with Lagrangian 
density £2 (PM1) is 


" = 1 1 ^ 

Let us define a new radial coordinate 


n • 2 , 

72 sm J 


-2 K s - sinV (| 


r dr dO 


(27) 


r = qr 


where q is a constant. Then 


d£_ 

dr 

and the mass per unit length (1271) can be rewritten as 


df 

= ^ q 


H = 



d£ 

dr 


+ — sin 2 / 


- 2 q 2 K s 30 sin 2 / 



r dr dO 


(28) 


(29) 


(30) 


From (1301) . it is clear that the energy per unit length g goes to —00 as q —» 00 . This is an indication that the 
vortex is unstable to uniform stretching, as it would be energetically favourable for the field profile / to dilate without 
bound and the vortex effectively to evaporate to infinity. A more rigorous statement of this property follows on from 
Derrick’s theorem^, which states that a necessary condition for vortex stability is that (d/i/dq)\ q -\ = 0. It is evident 
that (l30l) does not satisfy this criterion. 

In an attempt to fix this problem, we could add a ’’mass” term K v ( 1 — z ,cj>) to the Lagrangian density £ 2 , where z 
is the direction of 4> at r = 00 (where / = 0). The Lagrangian density 

£3 = £2 + K V {1 — n.fy (31) 

corresponds to the baby Skyrmion model in equation (2.2) ofii, p.207. The kinetic term together with the Skyrme 
term in £2 are not sufficient to stabilize a baby Skyrmion, as the kinetic term in 2 + 1 dimensions is conformally 
invariant and the baby Skyrmion can always reduce its energy by inflating indefinitely. This is in contrast to the 
usual Skyrme model, in which the Skyrme term prohibits the collapse of the 3 + 1 solitoni 2 . The mass term is added 
to limit the size of the baby Skyrmion. 
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IV. SKYRMION VORTEX WITH A TWIST: THE TWISTED SKYRME MODEL 

Instead of adding a mass term to stabilize the vortex, we will retain the Skyrme model Lagrangian (1241) but include 
a twist in the field g in (fill) . That is, instead of choosing 


g = nO-x 

we choose 

g = nO + mkz 

where mkz is the twist term. The Lagrangian density (1241) then becomes 


£2 = 


2A 2 


+ s i n / ( — r + 


2 / 2 
r k 


- 2 K s sin 2 / 


($_ 

\dr 


n~ 


2 ; 2 
i k 


(32) 

(33) 


(34) 


The value of the twist lies in the fact that in the far field, where r —> oo then / —> 0, the Euler-Lagrange equations for 
/ for both £3 without a twist and C 2 with a twist are formally identical to leading order, with m 2 k 2 / A 2 in the twisted 
case playing the role of the mass coupling constant K v . So it is expected that the twist term will act to stabilize the 
vortex just as the mass term does in £ 3 . 

On a physical level, the twist can be identified with a circular stress in the plane perpendicular to the vortex 
string (which can be imagined e.g. as a rod aligned with the 2 -axis). The direction of the twist can be clockwise 
or counter-clockwise. In view of the energy-mass relation, the energy embodied in the stress term contributes to the 
gravitational field of the string, with the net result that the trajectories of freely-moving test particles differ according 
to whether they are directed clockwise or counter-clockwise around the strings. 

The Euler-Lagrange equation (1251) corresponding to the twisted Skyrmion vortex Lagrangian density (I34[) reads 


0 = 


1 

A 2 


dr 2 
n 2 


1 # 

r dr 


■ m 2 k 2 


- 4 ( — + m k 


22 ' K s sin / cos / ( 


sin / cos / 

# x2 


-4 


— + m 2 k 2 ) K s sin 2 / 


dr 2 


1 df_ 
r dr 


(35) 


It should be noted that the second Euler-Lagrange equation ECU) is satisfied identically if g has the functional form 

(El. 


V. GRAVITATIONAL FIELD OF A TWISTED SKYRMION STRING 


We are interested in constructing the space-time generated by a twisted Skyrmion string. Without gravity, the 
Lagrangian density of the system is £ 2 , as given in equation El- 

To add gravity, we replace in C 2 with a space-time metric tensor g^ , which in view of the time-independence 
and cylindrical symmetry of the assumed vortex solution is taken to be a function of r alone. Metric tensor, g (lv , is 
of course the inverse of the covariant metric tensor, g MJ/ , of the space-time where g^ v = (ff Mt /) _1 - We use a cylindrical 
coordinate system (t, r, 9 , z), where t and 2 have unbounded range, r £ [0, 00 ) and 8 £ [0, 27r). The components of the 
metric tensor 


fgu 0 0 0 \ 

0 g rr 0 0 

0 0 g ge g Sz I 

V 0 0 g z6 g zz J 


(36) 


are all functions of r, and the presence of the off-diagonal components gg z = g z g reflects the twist in the space-time. 
The Lagrangian we will be using is 


U = 2^2 w d rf d »f + sin2 / 9^ d,,g d„g) 

+ 2 K s sin 2 / [(g^ v d^f d u f)(g KX d K g d x g) - 


with / = /(r) and g = nO + mkz. 
We need to solve 


2 sin 2 f ( 1 g^ d^f d„g) 2 } 


(37) 
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(i) the Einstein equations 

G - T 

where the stress-energy tensor of the vortex, Xis defined by 

rr - n dc * r 

1 Z Qgixu SW ^4 

and 


with R^ the Ricci tensor and 

the Ricci scalar; and 
(ii) the field equations for / and g 

and 


Gfj,v — Rfiv c^g^v R 


R = g„u R^ = g" 1 ' R„ 




v^ 1 


<9£ 4 

9£ 4 

d(df/daA) 

" df 

dU 

dC A 

d(dg/dx fl ) 

dg 


(38) 

(39) 

(40) 

(41) 

(42) 

(43) 


However, the field equations for / and g are in fact redundant, as they are satisfied identically whenever the Einstein 
equations are satisfied, by virtue of the Bianchi identities V^G() = 0. So only the Einstein equations will be considered 
in this section. 

To simplify the Einstein equations, we first choose a gauge condition that narrows down the form of the metric 
tensor. The gauge condition preferred here is that 

gee gzz - ( gez ) 2 = r 2 (44) 


The geometric significance of this choice is that the determinant of the 2-metric tensor projected onto the surfaces of 
constant t and z is r 2 , and so the area element on these surfaces is just r dr dd. 

As a further simplification, we write 


gtt — A ; g rr — 13 

gee = —C 2 ; ge z = w 


and so 


g zz = 


C 2 


The metric tensor, g^, therefore has the form 


9/j.v = 


(A 2 0 0 

0 —B 2 0 

0 0 -C 2 

\ 0 0 w 


0 

0 

U) 




(45) 

(46) 


(47) 


(48) 


Substituting (1351) into £4 gives 


£4 = {B~ 2 f' 2 + sin 2 /[n 2 (l + u> 2 /r 2 )C~ 2 + 2mkn r~ 2 uj + m 2 k 2 r~ 2 C 2 ]} 

+ 2 K s sin 2 / {R- 2 / ,2 [n 2 (l + uj 2 /r 2 )C~ 2 + 2mkn r~ 2 oj + m 2 fc 2 r“ 2 C 2 ]} 


( 49 ) 
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From equations (139)) . (l40l) . (l48l) and (1491) . the non-zero components of the stress-energy tensor, Tare: 

T tt = ^ A 2 B~ 2 f' 2 - A 2 + 2 K s B~ 2 f l2S j sin 2 f[n 2 { 1 + w 2 r“ 2 )G" 2 + 2 mkn wr“ 2 + (mfc) 2 C 2 r" 2 ] (50) 

T rr = JX 2 ^ “ B * ( 2^2 + 2K ° B ~ 2 f' 2 ^ sin2 / [ n2 (! + w 2 r _2 )C _2 + 2 mkn uor~ 2 + ( mk) 2 C 2 r ~ 2 ] (51) 

B~ 2 C 2 f' 2 + + 2 K s B~ 2 f' 2 ^ sin 2 / [n 2 (-l + w 2 r“ 2 ) + 2mfcra W C 2 F 2 + (mfc) 2 r- 2 G 4 ] (52) 


Bsz = ~ 2 X 2 B 2 w / /2 + (^^2 — 2K S B 2 /' 2 ^ sin 2 / [n 2 (l + w 2 r 2 )wC 2 + 2mkn(l + uj 2 r 2 ) + (mk) 2 r 2 wC 2 ] (53) 


and 


T = — 

1 rr 


2A 2 


( W 2 + r 2 )B~ 2 C~ 2 f' 2 + ( ~ + 2 K s B~ 2 f 2 ) sin 2 / 


[n 2 r 2 (l + w 2 r - 2 ) 2 G -4 + 2mfcn(l + r _ 2 w 2 )uC “ 2 + (mfc) 2 (-l + r _ 2 w 2 )] 


(54) 


where /' = df /dr. Note that Tg z is in general non-zero, provided that either mk or w is non-zero. In fact, mk acts as 
a source term for to, as will be seen in equation (l70l) below. The twist in the vortex is therefore solely responsible for 
a non-zero circular stress Tg z . 

In component form, the Einstein equations (13811 read 


Gtt — 

8nG 

G rr - = 

8ttG 

c 4 Ut ' 

T 

c 4 rr 

Gee = 

8ttG 

G* z = 

8nG 

,4 

,4 


(55) 

(56) 


and 


G zz = - 


87 tG 


T, 


(57) 


The non-zero components of the Einstein tensor, G^ Vl corresponding to the metric tensor (1401) are 


_ r) 1 A 2 B' A 2 C' A 2 C' 2 „ _ 2 o\ AWC' A 2 oj' 2 

G tt - Ra - ~gtt R - —-ttt ~ 2 ^ + t, 2^2 C 1 + r w ) ~ 


Gr 


n 1 A! C' G' 2 _2 2 , ww'C' 

- ^ rC + C 2^ 1+r W ) t 2 C 


//I 

4r 2 


(58) 

(59) 


Gee = 7?ee — R 
C 2 A! C 2 A" 


C 2 B' C 2 A'B' 2 CC CC" C' 2 , 2 Q\ 

+ - ^(2 + 3 r - 2 W 2 ) 


MB 2 AB 2 rB 3 AB 3 rB 2 
CA'C' CB'C' 3 ww'GG' 3GV 2 


AB 2 


I ? 3 


r 2 B 2 


4r 2 B 2 


(60) 


Ge z = i?e z — b 


ojA' 


UJ 


uA" 


ujB' w'B' ojA'B' 3 uC' 8 uj 2 u}'C' 


rAB 2 2 AB 2 AB 2 rB 3 2 B 3 AB 3 

V 2 ,.,// 3ww /2 


rB 2 G r 2 B 2 C 


3wC" _ 9 w ' 

52^2 (! + r w ) + 


w 


+ 


2 rB 2 2B 2 4r 2 B 2 


( 61 ) 

































with Gg z = G z g , and 


Gzz — Rzz i 2^ ZZ ^ 


Substituting 


oj 2 A' 


rAB 2 C 2 

oj 2 B' 

' r B^C 2 H 
uu'A! 


r 2 A" 
AB 2 C 2 
r 2 B'C' 


r 2 A' n' r 2 A!(-if 

(-1 I — 2 2\ , ' ^ JD (-1 , -2 2\ l ^ — 2 2\ 

(! + r ^ „ 0.^2 ( 1 + r u )- u,^ ( 1 + r w ) 


B 3 C 3 

ww'S' 


(1 + r 2 w 2 ) + 


AB 3 C 2 
4 w 2 C" 


■ r w 
wu'C' 


rB 2 C 3 
(1 - 3r" 2 w 2 ) - 


B 2 C 3 

LULU 1 


AB 2 C 3 

o, .2(~il2 

(l + r - 2 w 2)_^_( 1 + r -2 fa; 2 ) 




B 2 C 4 

/2 


W 


AB 2 C 2 B 3 C 2 B 2 C 3K ' rB 2 C 2 

(l62l) and (l50l) - (l54l) into equations (f55]) - ((57]) now gives 

A 2 C' A 2 C' 2 


B 2 C 2 AB 2 C 2 


(1-3 r~ 2 u 2 ) 


8t tG A 2 B' 

—Ttt — - 


rB 3 rB 2 C B 2 C 2 


(1 + r~ 2 uj 2 ) — 


A 2 


'C' 


8t tG 


rr rA rC C l{ 


A 2 oj' 2 


r 2 B 2 C Ar 2 B 2 


f 2 


_<> .. x V_X UJ 

+ r CU ) - ry -h —2 

r 2 C 4r 2 


( 62 ) 

(63) 

(64) 


87tG 


Tgg = — 


+ 


C 2 ^' C 2 A n 


rAB 2 
CA'C' 


AB 2 rB 3 


C 2 B' C 2 A'B' 2 CC' 

+ -7T75— + 


AB 3 


rB 2 


/-i/-iff /-if 2 

UU vs . _o 9 

_ B 2__ S^ (2 + 3r W 


8t tG 


Tg z = 


AB 2 

B 3 

r 2 B 2 Ar 2 B 2 




uiA' 

w'A' 

luA" 

wB' u/B' 

loA’B' 

3 wC" 

3u> 2 uj'C' 

rAB 2 

2 AB 2 

AB 2 

rB 3 2 B 3 

AB 3 

rB 2 C 

r 2 B 2 C 

3 loC' 2 

(1 + r 

-V) + 

w' w" 

3wu/ 2 



B 2 C 2 

2rB 2 2B 2 ' 

4r 2 B 2 




(65) 


( 66 ) 


and 


8t tG 


u> 2 A' 


„ r 2 4" 

Tzz = ~ rAB 2 C 2 AB 2 C 2 


r~A" _ 2 2 r 2 A'B' _ 2 2 r 2 A'C" _ 2 2 

(! + r ^ ) + -n5s7«( 1 + ^ ^ )-T5527rl( 1 + ^ ^ ) 


AB 3 C 2 


to 2 B' r- 

+ 


-(1 + r~ 2 w 2 ) + 


rB 3 C 2 B 3 C 3 
ojlo'A' ww'B' ojlu'C 1 


rB 2 C 3 B 2 C 3 


AB 2 C 3 

o. ,2/~u2 

(1 + r -W)-^-(l+r-W) 


AB 2 C 2 B 3 C 2 B 2 C 3 


(1 - 3 r _ 2 or) - 


ww 


+ 


B 2 C 4 

12 


rB 2 C 2 B 2 C 2 4B 2 C 2 


(1 - 3 r _ 2 ar) ( 67 ) 


These equations can be rearranged as source equations for the second derivatives of A, C and ui and the first 
derivatives of B and / as follows: 


A' 


A" = - — 


A'B ' AB' AC' ojuj'AC' AC ' 2 


+ 


Au> 


12 


C" = - 


2 r ' B ’ 2rB ' rC ' r 2 C 
C' CA’ A'C' CB' B'C' 


2 rA 


2 rB 


B 




2 UJUj'C' C' 2 „ 2 2, Cuj' 2 

__ + _ {1 + 2r - W) + __ 


- KB 2 C- l r~ 2 [(nu + C 2 km) 2 - (nr) 2 ] ( + 2K S B~ 2 f' 2 ) sin 2 / 


1 


( 68 ) 


(69) 


UJ 

r 


ujlj' 2 ujA! 

r 2 rA 


4 uC' 4 wVC" 4wC" : 


A 


rB 


B 


rC 


r 2 C 


C 2 


(1 + r _2 w 2 ) 


-2 KB 2 C~ 2 r - 2 {u[(noj + C 2 km ) 2 + (nr) 2 ] + 2C 2 kmn r 2 } + 2K s B- 2 f 2 ^j sin 2 / 


(70) 


S' = A~ 1 BA' + -^B 3 C _2 r _1 [(nu + C 2 fcm) 2 + (nr) 2 ] sin 2 / 


(71) 
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and 


r = 


-l 


a' a c ' 2 

Va + Vc~c^ { 


2, ,2 


) + 


r 2 C 


4 r 2 


+ -^S 2 C- 2 r- 2 [(na; + C 2 fcm ) 2 + (n«) 2 ] sin 2 / 


/— 

\ 2A 2 


— 2K S C 2 r 2 [(nw + C 2 km) 2 + (nr) 2 } sin 2 / 


(72) 


where k = 87 rG/c 4 . 

In order to solve equations ©-ED, we require boundary conditions on the functions >1, 13, G, oj and / for small 
r and in the limit as r tends to 00 . In particular, ui should vanish both at r = 0 (at least as rapidly as r 2 , so as to 
preserve elementary flatness on the axis of symmetry) and in the limit as r —> 00 (to allow for a locally Minkowski 
vacuum at space-like infinity). It is evident from equation (1701) that u> will be zero everywhere unless mk 7 ^ 0, and so 
the vortex twist acts as a source for the space-time twist as mentioned above. 


VI. THE EINSTEIN-HILBERT ACTION 

The Lagrangian density £4 of course describes the matter part of the twisted Skyrme model only. When we vary 
£4 with respect to / and g, we obtain the Euler-Lagrange equations of motion for the twisted Skyrme model, and 
when we vary £4 with respect to g^ we obtain the stress-energy tensor for the twisted Skyrme model. But these are 
not in themselves sufficient to generate the gravitational field equations for g AtJ/ . 

To do this we need to add to £4 the Hilbert Lagrangian density, which is proportional to \J (—A )R, where A is the 
determinant of the metric tensor, g^, and R is the Ricci scalar. In Riemannian geometry, the Ricci scalar or scalar 
curvature is the simplest curvature invariant on a Riemannian manifold. It assigns to each point on a Riemannian 
manifold a single real number determined by the intrinsic geometry of the manifold near that point. Specifically, 
the scalar curvature represents the amount by which the volume of a geodesic ball in a curved Riemannian manifold 
deviates from that of the standard ball in Euclidean space. In two dimensions, the scalar curvature is twice the 
Gaussian curvature, and completely characterizes the curvature of a surface. 

In general relativity, the Ricci scalar acts as the Lagrangian density for the Hilbert action. The Euler-Lagrange 
equations for this Lagrangian density under variations in the metric tensor constitute the vacuum Einstein field 
equations. In other words, if £4 = 0 (or if k = 0 so that there is no coupling between gravity and matter) then the 
Hilbert action will give us the vacuum Einstein field equations. 

The Einstein-Hilbert action 

I = J \J ( — g)(R + k £4) d 4 x ( 73 ) 

£5 

(or equivalently the Lagrangian density £5) is a functional of two sets of fields: the Skyrmionic string fields / and g, 
and the metric tensor components g^. But only £4 depends on / and g explicitly. This means that £4 depends on / 
and g , but R does not. Of course both £4 and R depend on g lliy . So when we vary £4 with respect to / and g , we can 
forget R (and yj— A) and just vary £4 with respect to / and g. This gives us the equation of motion of the string. 

When we want to calculate the effect of the stress-energy content of the string on the gravitational field around it, 
we need to vary the Einstein-Hilbert action with respect to g or more precisely the inverse metric tensor g ^ v . This 
means that we need to vary all the terms in £5 with respect to g ^ v , including R, £4 and \/—A. This then gives us 
the Einstein equations, which are second-order PDEs for the components of g^. 
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